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Evidence for a relativistic iron line in GRS 1915+105 
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Abstract. We report on the discovery of a relativistic iron Ka fluorescent emission line in the BeppoSAX spectrum 
of the microquasar GRS 1915+105, taken on April 19, 1998, when the line was unusually intense. The feature 
is broad and skewed, clearly indicating emission from the innermost regions of the accretion disc. The inner 
emitting orbit is larger than the innermost stable orbit, even in Schwarzschild metric: thus, a non-zero BH spin 
is not required by these data. 
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1. Introduction 

1.1. Iron line diagnostics of relativistic accretion discs 

Black Hole (BH) accreting systems are expected to exhibit 
spectral features whose form is determined by reprocess- 
ing of primary high-energy photons by gas swirling around 
the central body. Resulting spectral distortions are of ut- 
most importance because they probe the effects of strong 
gravity, which are otherwise difficult to reveal. General- 
relativity (GR) calculations of predicted spectra, and their 
comparison with observational data, may help to mea- 
sure physical parameters of BHs, such as their masses and 
angular momenta (e.g. Martocchia 2000, and references 
therein). 

The iron Ka fluorescent emission is specially impor- 
tant as a diagnostic of GR effects (e.g. Fabian et al. 1989; 
Laor 1991; Matt et al. 1992; Martocchia & Matt 1996; 
Martocchia et al. 2000). In case it is originated in the 
innermost regions of an accretion disc, it gives the pos- 
sibility, in principle, of distinguishing between a rotat- 
ing (Kerr) and a static (Schwarzschild) central BH. This 
is because the resulting line profile is affected both by 
Special Relativity (Doppler, including transverse, shifts) 
and General Relativity (energy shift, lcnsing) effects. It is 
very broad and asymmetric, the exact shape critically de- 
pending on the geometrical parameters of the emitting re- 
gion (disc inner and outer radii, inclination angle). For the 
same set of parameters, differences between Schwarzschild 
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and extreme Kerr line profiles are rather subtle. However, 
the innermost stable orbits of the disc are different in the 
two cases - corresponding respectively to 6 and 1.23 grav- 
itational radii - and the profiles are very different for the 
two values. The measurement of an innermost emitting 
radius less than 6r g can therefore be the signature of a 
spinning BH. 

This diagnostic has been successfully used in the case 
of Seyfert galaxies (e.g. Tanaka et al., 1995, Wilms et al., 
2001), while in galactic BH candidates ionization of the 
inner accretion flow generally prevents fluorescent emis- 
sion from being effective. Nevertheless, intense, broad Fe 
Ka emission has been observed in Cyg X-l with both 
BeppoSAX (Frontera et al., 2001) and Chandra (Miller 
et al., 2002a). Large, possibly relativistic Ka profiles have 
been detected also in XTE J1748-288 (Miller et al., 2001) 
and XTE J1650-500 (Miller et al., 2002b). 

1.2. GRS 1915+105 

GRS 1915+105 is one of the so-called microquasars, 
Galactic jet sources with properties similar to those of 
quasars, but on a smaller, stellar scale. It was extensively 
studied since its discovery by the WATCH experiment on 
GRAN A Tin 1992 (Castro-Tirado et al. 1992). Due to very 
large interstellar absorption, the optical counterpart of 
this source was not easily found: only recently, Greiner et 
al. (2001a) reported results of a spectroscopic analysis in 
the H and K bands, which suggest that the mass-donating 
star is a K-M III star, i.e. that GRS 1915+105 belongs to 
the class of low-mass X-ray binaries. 
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Greiner et al. (2001b) used these infrared spectro- 
scopic observations to infer the companion's orbital period 
- which is ~ 33.5 days - and determine the mass function 
f(M c ) ~ 9.5Mq. Since the donor's mass, for K-M III stars 
of that kind, is known to be ~ 1.2Mq, and the inclination 
is known, too, the mass of the central compact object has 
been constrained to M c = 14 ± 4M Q , i.e. well above the 
standard neutron star mass limit. 

Weak variable radio/IR emission connected with jets 
has been often detected (Mirabel & Rodriguez 1994). 
Measurements of the jets' superluminal motion allowed 
to estimate the source distance D ~ 12.5 kpc and incli- 
nation i ~ 70° (Rodriguez & Mirabel, 1999; Fender et al., 
1999). 

GRS 1915+105 is well-known for the extremely wide 
variety of its variability modes, ranging over all possible 
bands and timescales (e.g. Markwardt et al. 1999, Belloni 
et al. 1997a, b; Greiner et al. 1996). Time- resolved obser- 
vations have been intensively performed especially with 
the Rossi X-ray Timing Explorer (RXTE) satellite mis- 
sion. Belloni et al. (2000), for instance, managed to clas- 
sify these observations into 12 separate classes based on 
count rate and color characteristics. In this way, the source 
variability could be explained in terms of transitions be- 
tween three basic states: two soft states with a fully ob- 
servable disc at two different temperatures ( "A" with the 
lower temperature, "B" with the higher), and a hard state 
corresponding to the non-observability of the innermost 
parts of the accretion disc (state "C"). State "C", with 
low flux and hard color, can last uninterruptedly for weeks 
or even months. The most popular model for the source's 
large scale variability is the one proposed by Belloni et al. 
(1997a, b), who invoked the onset of a Lightman-Eardley 
thermal-viscous instability which blows-off the innermost, 
radiation-pressure supported part of the (optically thick) 
disc. The hard state (C) would then correspond to a steady 
replenishment, which takes place on a viscous timescale. 
Indeed, evidence of the temporary disappearance and sub- 
sequent restoring of the inner portion of the accretion disc 
comes from simultaneous BeppoSAX and Ryle radiotele- 
scope observations (Feroci et al., 1999). 

The iron Ka emission line in GRS 1915+105 is usu- 
ally not very intense, and is therefore difficult to study in 
detail - not to mention the dramatic spectral variability 
of the source, which requires selecting short time inter- 
vals, so reducing the statistical quality of the spectra. Lee 
et al. (2001) recently reported that neutral Fe Ka emis- 
sion is likely present in the high-resolution data of a low- 
hard state, collected with Chandra HETGS, together with 
Fe photoelectric edges and absorption lines due to highly 
ionized matter. These features indicate that ionized iron, 
probably belonging to the accretion disc itself, is present 
on the line of sight with abundance larger than solar. 

Actually, during the several observational campaigns 
performed by BeppoSAX, the Fe Ka line was almost in- 
variably rather weak, with Equivalent Widths (EW) of a 
few tens of eV at most. Only in the observation of April 
19, 1998, the line was strong enough (EW > 200 eV; cp. 
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Fig. 1. The MECS lightcurve of GRS 1915+105 during 
the April 19, 1998 BeppoSAX observation, which lasted ~ 
120 ks. The time intervals we considered are those marked 
(comprised between t ~ 3 x 10 4 and t ~ 6.2 x 10 4 s). In 
this paper, the fifth of these intervals is discussed in more 
detail. 



Fig. ||) for a detailed analysis of its profile to be performed. 
We have therefore analyzed the spectral data from this 
observation, further selecting some intervals of relative 
calmness in the lightcurve (see Fig. [l]), to limit the spec- 
tral variability as much as possible. In order to model the 
source emission, we adopted the illumination model, local 
spectra and emissivities described in Martocchia, Matt & 
Karas (2001) and Martocchia, Karas & Matt (2000; see 
also: Martocchia 2000). Finally, we used our own fitting 
tool, dedicated to GR spectral distortions in Kerr met- 
ric, to analyze the data. The tool, named kerrspec, is 
described in the Appendix. 

2. Data reduction 

BeppoSAX (Boella et al 1997) is composed by several in- 
struments. Here we made use of the two imaging instru- 
ments, the LECS and the MECS, and the high energy, 
passively collimated instrument, the PDS. As customary, 
we have used the following energy bands: 0.1-4 keV for 
the LECS; 1.8-10 keV for the MECS; 15-150 keV for the 
PDS (neglecting higher energy data because of insufficient 
statistics due to the short integration times). The LECS 
and MECS spectra have been selected from regions of 8' 
radius and centered on the source. The background has 
been evaluated from similar regions in blank field obser- 
vations, and then subtracted. The PDS is composed by 
four units, two of which (switching every 96 s) monitoring 
the background, which is automatically subtracted in the 
pipeline producing the spectrum. 

We discuss here the BeppoSAX observation of 
GRS 1915+105 started on April 19, 1998 at 11:33:19 UT 
and ended on April 20, 1998 at 20:14:52 UT. Among the 
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several BeppoSAX campaigns on GRS 1915+105, we se- 
lected this one because the iron line EW was significantly 
larger than in all other observations. 

Out of the entire ~ 120 ks observation, some time in- 
tervals were sorted out where the source variability is less 
pronounced. These are the six satellite orbits starting from 
3xl0 4 s after the beginning of the observation, see Fig. 1. 
Analysis of the other parts of the observation would re- 
quire a finer temporal analysis which is beyond the scope 
of this paper and is deferred to a future work. 

All spectral fits were performed with the xspec v. 11 
software package. In the following, all errors refer to 90% 
confidence level for one interesting parameter (A% 2 =2.71) 

At the flux level of GRS 1915+105, it is very likely 
that systematic errors cannot be neglected with respect 
to the statistical ones. This may explain the rather large 
values of the \ 2 reported below. 

3. Data analysis and results 

3.1. The best-fit model in a selected time interval 

Let us first discuss the fifth interval as an example, being 
the one in which a low value of \ 2 was found when fitting 
with simple models. The model for the continuum includes 
the usual ingredients for this source: a disc thermal emis- 
sion (diskbb), a powerlaw with exponential cutoff, and 
a Compton Reflection component (pexrav: Magdziarz 
& Zdziarski 1995), which is expected to always accom- 
pany the iron emission line. All these components are seen 
through interstellar absorption. The best-fit spectrum and 
A\ 2 for the fifth interval, fitting without any iron line, is 
shown in Fig. |[ left panel (x 2 cd = 2.64). The presence 
of an iron line is obvious. To better show the profile of 
the line, we also fitted the spectrum excluding the energy 
region which contains this feature (3-7.5 keV), and then 
reintroduced this energy interval. The residuals are shown 
in Fig. H (fifth panel), after reinserting the data of the 3- 
7.5 keV interval. The presence of a prominent, broad and 
red-tailed iron line is clear. 

We first tried to model the line with a gaussian. The 
best-fit is obtained with a rest energy Eq — 6.38 ± 0.12 
keV, a = 0.70 ± 0.15 keV and EW- 250 eV; the x? cd is 
then 1.48 (see Fig. ||, central panel). The most natural 
interpretation for a so broad line is emission from a rel- 
ativistic accretion disc (Fabian et al., 1989), an hypoth- 
esis further supported by the asymmetry of the profile 
(see Fig. ^ . The next step was therefore to substitute the 
gaussian with the relativistic model kerrspec, described 
in the Appendix. To be self-consistent, we also applied the 
relativistic corrections to the reflection continuum (using 
kerrspec) as well as to the disc thermal component (us- 
ing diskpn, Gierlinski et al. 1999). After having verified, 
however, that the results are insensitive to the latter cor- 
rections, in order to save computer time we adopted the 
simple, non relativistic models for the continuum in all the 
following fits. The inclination angle has been always fixed 
to 70° (Mirabel & Rodriguez 1994), while the inner and 



outer disc radii were kept free to vary. We tried different 
emissivity laws, and found the best solution with a rather 
"flat" dependence on radius (oc r~ 2 ). With the relativis- 
tic line model (Eq = 6.4 keV) the fit is better than with 
a gaussian: Xr<x\ = 1.31. The spectrum and residuals are 
shown in Fig. 0, right panel. If Eq is instead fixed to 6.7 
keV, corresponding to He-like iron, the fit is significantly 
worse (x 2 cd » 2). 

We also tried adding a narrow gaussian centered either 
at 6.4 or at 6.7 keV; no significant improvements in the 
fits were found. 

One problem with the above fits is that the resulting 
parameters of the thermal emission are rather unplausi- 
ble. The temperature is very large, kT = 3. 77^43 keV, 
and the component's normalization very low (0.41 ± 0.22; 
see Arnaud 1996 for the meaning of this normalization), 
implying a value of the inner radius of the order of one 
kilometer. This, in turn, implies an upper limit to the BH 
mass (obtained assuming the inner emitting radius con- 
doling with the innermost stable orbit for a maximally 
rotating BH) of less than 1 solar mass, contrary to ob- 
servations (Greiner et al., 2001b). If the disc blackbody 
component is neglected in the model, the fit e.g. with the 
gaussian line yields a worst fit Xied = 1-61, but the re- 
sults on the line parameters does not change significantly. 
Moreover, the fit with a relativistic line is still better than 
with a gaussian also in this case. 

We therefore conclude that the relativistic solution for 
the iron line, besides being the most natural physical ex- 
planation for the very broad line observed, does not de- 
pend on the details of the continuum model. 

The small relative luminosity of the disc thermal com- 
ponent (if present at all) brings us to interpret the source's 
state as a 'C'-type, in the scheme of Belloni et al. (2000). 
This interpretation is in agreement with the occurrence of 
a "plateau" in the light curve (lower flux and variability), 
although a hard spectrum would be expected, while in- 
stead we get a steep powerlaw, and a cutoff at rather low 
energies (see Table 1.). 

3.2. Other time intervals 

For all other time intervals we found qualitatively similar 
results: a broad and prominent iron line, which is better 
fitted by a relativistic line profile than by a gaussian, and 
a faint disc blackbody component. 

An inspection of the residuals often indicates prob- 
lems both around 1 keV, and above 50 keV. The latter 
problem may be due to a non-thermal component aris- 
ing when the Comptonized power law cuts-off, while we 
have no obvious explanation for the former problem. To 
overcome these difficulties, and because we are interested 
here to the iron line properties, we decided from now on 
to exclude the LECS data in the fits, and to use the PDS 
up to 50 keV only. With this prescriptions, and assum- 
ing as working hypothesis that the iron line arises from a 
relativistic disc, we analysed all the time intervals with a 
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Fig. 2. Left: The best-fit spectrum and A\ 2 for the fifth interval, fitting without modelling the iron line (x 2 c d = 2.64): 
the residuals are evident. Center: The best- fit spectrum and A\ 2 for the same interval, when the line is modelled by 
a gaussian profile (x 2 e d = 1-48). Right: The best- fit spectrum and A% 2 for the fifth interval again, with a fully-GR 
model of the line profile (x 2 ed = 1-31). See the text for details. 




Fig. 3. Residuals at the iron line energies for all the selected time intervals. For illustration purposes, the fits have 
been performed excluding the energy region of the iron line. Prominent, broad and skewed features are then evident 
around ~ 6.4 keV. 



model which neglects disc thermal emission and includes, 
besides the line, a power law with exponential cut-off and 
a Compton reflection continuum. We stress again that the 
latter component is included more on astrophysical ground 
(it is supposed to always go along with the iron line) rather 
than because strongly required by the data. After having 
verified that the values of the absorbing column densities 
riH for the various intervals are all consistent with each 
another and to the value of 5.4xl0 22 cm -2 , for simplicity 
we fixed in all fits nu to this value. 

The results for the various time intervals are summa- 
rized in Table [j] The larger value of the x 2 e d 01 tne nltn 
interval with respect to those described above is due to 
the lower number of spectral bins and to the exclusion of 
the thermal component. Given the abovementioned un- 
certainties on the proper modeling of the continuum, we 
stress that the values in the table should be taken with 
caution: a different continuum models may result in some- 
what different values of the disc parameters. In any case, 



a general result seems to emerge: the inner radius of the 
line emitting region is larger than the innermost stable 
orbit. Whether this means that the innermost regions are 
absent, or simply that they do not emit the iron line (for 
instance because the matter is mildly ionized and resonant 
trapping effective, Matt et al. 1993) we cannot say. 

Finally, we tried summing the spectra of all six time 
intervals together, and fitting with the same model. The 
very large reduced \ 2 , 5.1, indicates that spectral variabil- 
ity is too large over these time scales. In particular, strong 
residuals at the iron line energies imply that the line pa- 
rameters change with time, as suggested by the results of 
the time-dependent spectral analysis. 

4. Discussion 

A strong iron Ka fluorescent line is present in the April 19, 
1998, spectrum of the superluminal galactic source GRS 
1915+105. The line is broad and skewed, clearly indicat- 



A. Martocchia et al.: Evidence for a relativistic iron line in GRS 1915+105 



5 



Table 1. The best- fit parameters for the various time intervals we considered, excluding LECS data and using a 
simplified model, which does not take into account thermal emission from the disc (wabs ( PEXRAV + kerrspec ) ) . 
The rest line energy is 6.4 keV. At is the MECS exposure time. See the text for more details. 
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ing emission from a relativistic accretion disc. We success- 
fully fitted the line with a fully relativistic profile. To our 
knowledge, this is the first time that a relativistic iron line 
is observed in this source, largely due to the fact that we 
catched the line with an unusually large equivalent width. 

While providing evidence that an accretion disc resides 
in the innermost part of this system, our result confirms 
that strong gravity is at work in distorting the iron line 
profile. From the best-fit parameters we conclude that 
the line is from neutral or low-ionization iron, and likely 
emitted from a region of the accretion disc which is close 
but probably does not include the innermost stable orbit, 
even for a static BH (but the values of the disc parameters 
are somewhat dependent on the adopted model for the 
continuum, and must be taken with caution). This may 
imply either that the disc does not extend down to this 
radius, or that the line is not emitted in the innermost 
regions, maybe because the iron there is in the intervals of 
ionizations in which resonant trapping efficiently destroys 
line emission (e.g. Matt et al. 1993; 1996). Unfortunately, 
at these radii differences in the line profile between Kerr 
and Schwarzschild metrics are too subtle to be detectable, 
and the question of the angular momentum of the BH - 
which, from this analysis, is not required to have a large 
value - actually remains open. 

So, no final statement about the spin of the BH in 
GRS 1915+105 can still be made. Zhang et al. (1997), 
on the basis of soft state X-ray brightness (which is most 
likely due to thermal emission), suggested that a fast ro- 
tating BH is hosted at the center of this source. Their con- 
clusion was based on a mass value twice that measured 
by Greiner et al. (2001b). However, even accounting for 
the correct mass, the BH spin seems to be larger in GRS 
1915+105 than in GRO J1655-40, and this does not fit to 



the results based on QPO diagnostics, as explained below. 
The accretion disc temperatures of the two sources are 
nearly identical, and the mass estimate for GRO J 1655- 
40 is M ~ 7M Q (Orosz k Bailyn, 1997; Shahbaz et al., 
1999). Taking Zhang et al.'s innermost thermally emitting 
orbit for GRS 1915+105 (at ~ 40 km, i.e. r ms - 2r g ), one 
infers a BH spin a ~ 0.95 in this source (Bardeen et al., 
1972). Actually, Zhang et al.'s estimates of the BH angu- 
lar momentum are based on rather simplified models of 
the continuum, and they do not include GR spectral dis- 
tortions and effects of radiation transfer through the disc 
atmosphere (e.g. Merloni et al., 2000; Gicrlinski et al., 
1999 and 2001). In our data analysis, even accounting for 
GR effects on the thermal component, we got disc temper- 
atures and luminosities which are physically unplausible, 
and do not allow us to test assumptions regarding either 
the BH mass nor its spin. Non-standard disc models (e.g. 
slim-disc) could actually reproduce high temperature discs 
for non-rotating BHs. 

QPO diagnostics could therefore be useful as indepen- 
dent tests. Various models try to explain the occurrence 
of QPOs, suggesting that the oscillations correspond to 
the motion on the last stable orbit (Shapiro & Teukolsky, 
1983), and sometimes invoking GR effects such as "frame 
dragging" and relativistic precession (e.g. Merloni et al., 
1999; Stella et al., 1999). An interesting possibility is that 
the highest QPO frequency changes in connection with 
the variation of the disc innermost radius. 

QPOs in GRS 1915+105 have been observed in a wide 
range of frequencies (0.001 + 67 Hz, e.g. Morgan et al. 
1997). In GRS 1915+105 the QPO highest frequency (67 
Hz), unlike all the others, does not change much with time, 
and may therefore be related to the accretion disc inner- 
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most stable orbit r ms (Morgan et al. 1997). A similar case 
is given by the 300 Hz QPO in GRO J1655-40. 

Recently, a second high-frequency QPO has been de- 
tected in GRS 1915+105 (Strohmayer 2001), with v ~ 40 
Hz. In terms of the frame-dragging model, this lower fre- 
quency peak could be associated to the radial epicyclic 
frequency f ra d through the formula vq = v^cp — ^rad 
(the nodal precession can be connected to a further low- 
frequency peak, already detected at ~ 0.9 Hz by Morgan 
et al., 1997). Strohmayer also reports that the 67 Hz QPO 
did in fact drift in frequency by as much as 5 percent. If 
this frequency is associated with the Keplerian motion at 
the last stable orbit, then - assuming the BH spin is the 
same in GRS 1915+105 and GRO J1655-40 - it should 
scale with M c _1 . Thus, a factor ~ 2 + 3 of discrepancy 
would remain with respect to Greiner et al.'s estimate. 
This indicates that the BH spin in the two sources is dif- 
ferent. If the QPOs in GRS 1915+105 and GRO J1655-40 
are seen in the framework of the frame dragging scheme 
(e.g. Cui et al., 1998, for details on the model), and using 
Greiner et al.'s mass estimate, we have a <~ 0.8 and 0.95 
for the two sources, respectively. On the other hand, if 
the QPO frequency had a diskoseismic origin (e.g. Nowak 
et al., 1997), the model would require a maximally spin- 
ning BH in GRO J1655-40, and an almost static BH in 
GRS 1915+105 (cp. Greiner et al., 2001b, and references 
therein). 

Clearly, further tests for the BH spin would be of ex- 
treme importance. Iron line diagnostics may be crucial to 
solve the puzzle: better resolution and sensitivity spectral 
data, provided by satellites of the new (e.g. XMM-Newton) 
and future (e.g. Constellation- X, XEUS) generations, will 
be probably decisive in this respect. 
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Appendix: The computational technique 

We have computed GR effects acting on both the iron 
Ka line emission and the Compton reflection continuum. 
We adopt the usual assumption of a geometrically thin, 
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optically thick, cold (with respect to X-rays reprocessing) 
reflector extending over the equatorial plane. We take into 
account the directional anisotropy of the illuminating flux, 
caused by aberration, Doppler effect and light focusing. 
GR anisotropy of the illuminating flux is crucial to dis- 
tinguish static from spinning BHs, because the primary 
illumination, as well as the induced emissivity of the reflec- 
tion component, carry information about the spacetime in 
the very vicinity of the BH. The iron line profile and, in 
particular, its intensity and EW are affected by spacetime 
distortions (Martocchia & Matt, 1996; Martocchia et al., 
2000). 

The observed line profiles are often obtained from the 
data by subtracting the Compton-reflected component 
without taking into account important relativistic effects, 
such as the iron edge smearing, which may be relevant if 
the emission occurs down to the innermost stable orbit, 
in particular for spinning BH. In our computations (e.g. 
Martocchia et al., 2000), besides using emissivity laws de- 
rived from a proper GR treatment of the illuminating flux, 
we treated at the same time the GR distortions on the 
iron line together with the underlying continuum hump 
(Martocchia et al., 2000; Martocchia et al., 2001). 

In order to fit real spectra, a ray-tracing code was 
developed in the form of a user-defined subroutine, 
kerrspec (Martocchia 2000), which has been linked to 
the XSPEC data analysis package (Arnaud 1996). In the 
adopted scheme, a large set of light rays (null geodesies) 
in the Kerr BH spacetime is pre-calculated and stored. 
This allows an efficient access during the fitting loop. The 
code allows the radial, angular and energy dependencies 
of the local emissivity to be specified on input in a tab- 
ular or semi-analytical form. After fixing the desired grid 
resolution in the image plane, the required light rays cor- 
responding to all image pixels are retrieved. 

The original approach (Karas et al. 1992) was refined 
in order to improve the resolution around the inner edge of 
the disc, very near to the BH horizon, and to ensure bet- 
ter confidence in the region where strong focusing occurs. 
This turned out to be necessary because the computed line 
profiles are very sensitive to the radius of the innermost 
emitting orbit, especially if the emissivity varies rapidly 
with radius. We recall that the innermost emission area 
is usually identified with a ring located at the marginally 
stable orbit, ri«r ms , the value of which depending on the 
BH angular momentum. 

The local emissivity law, F(r), must be specified before 
running the code. Different choices are available: a power- 
law dependence F q (r) oc r~ q ; a "standard" Novikov- 
Thorne-Page model (Novikov et al. 1974); or the effect 
of illumination from a central source (Martocchia & Matt 
1996). In the latter case, h (the primary source height) 
is a phenomenological parameter which determines the 
form of the local emissivity law: small h corresponds to en- 
hanced concentration of primary X-rays near the horizon, 
resulting in substantial anisotropy of the reflected com- 
ponent. The local emissivity has an angular dependence, 
too, including possibly a limb-darkening law. This may 



be expressed either by the functional forms for fluores- 
cent and Compton-reflected emission given by Ghisellini 
et al. (1994), or, alternatively, by the numerically-derived 
tables from the Montccarlo computations by Matt et al. 
(1991), which yielded the whole spectra of the reflected 
component. 

The local profiles of the emissivity were interpolated by 
polynomials to save time machine (cp. Martocchia 2000, 
Martocchia et al., 2001). These emissivity laws are valid 
between the inner disc edge (assumed to be equal to or 
greater than the innermost stable orbit) and the outer 
edge. Because the emission tends to be concentrated in 
the central parts of the disc, the inner emitting radius is 
usually better constrained in the fits than the outer radius. 

Typically, the code may take up to a few hours, on a 
common DECalpha workstation, to fit a spectrum using 
curved geodesies in the very vicinity of the BH, i.e. ac- 
counting for radiation emitted at distances of the order of 
r e from the event horizon. 



